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Orientation of x-lines in asymmetric magnetic reconnection — 
mass ratio dependency 

Yi-Hsin Liu,^ M. Hesse,^ and M. Kuznetsova^ 

Abstract. Using fully kinetic simulations, we study the x-line orientation of magnetic 
reconnection in an asymmetric configuration. A spatially localized perturbation is em¬ 
ployed to induce a single x-line, that has sufficient freedom to choose its orientation in 
three-dimensional systems. The effect of ion to electron mass ratio is investigated, and 
the X-line appears to bisect the magnetic shear angle across the current sheet in the large 
mass ratio limit. The orientation can generally be deduced by scanning through corre¬ 
sponding 2D simulations to find the reconnection plane that maximizes the peak recon¬ 
nection electric field. The deviation from the bisection angle in the lower mass ratio limit 
can be explained by the physics of tearing instability. 


1. Introduction 

Magnetic reconnection is arguably one of the most im¬ 
portant energy conversion and plasma transport processes 
in solar and space plasmas. Among other effects, it deter¬ 
mines the energy entry from the solar wind into Earth’s mag¬ 
netosphere, and it enables energy transport and dissipation 
therein [Dungey, 1961]. At Earth’s magnetopause, reconnec¬ 
tion proceeds asymmetrically between magnetosheath plas¬ 
mas, namely solar wind plasmas compressed by Earth’s bow 
shock, and magnetospheric plasmas. The magnetosheath 
side has a typical magnetic field strength ~ 20 nT, density 
~ 15 cm“^ and plasma-~ 2; The magnetosphere side has 
magnetic field strength ~ 60 nT, density ~ 0.5 cm“^ and 
plasma-yd ~ 0.1 (e.g., Phan and Paschmann [1996]). The 
magnetic field shear can be an arbitrary angle 0. Consider¬ 
ing a planar current sheet, the x-line could develop at any 
angle from 0 to 0, where the fields in the plane normal to 
this orientation have opposite signs, as suggested by Cowley 
[1976]. It is unclear if there is a simple principle to de¬ 
termine the orientation of the x-line in a three-dimensional 
(3D) system. 

The first attempt to address this fundamental problem 
was by Sonnerup [1974] (also independently by Gonzalez and 
Mozer [1974]), who suggested that reconnection will occur 
in a plane where the guide field is uniform. Motivated by 
Cowley [1976], the angle that bisects the total shear has been 
employed in global modeling [Moore et al., 2002; Borovsky^ 
2008; Sibeek^ 2009]. More recently, other sophisticated mod¬ 
els based on maximizing various physics quantities were pro¬ 
posed. Swisdak and Drake [2007] suggested the plane in 
which the reconnection outflow jets have a maximum speed. 
Sehreier et al. [2010] pointed out another possibility by max¬ 
imizing the reconnection electric field (equivalent to the re¬ 
connection rate), where the formulation in Cassak and Shay 
[2007] for asymmetric reconnection could be used. Based 
on 2D simulations at different oblique reconnection planes, 
Hesse et al. [2013] further proposed that the x-line orienta¬ 
tion should be determined by maximizing the peak recon¬ 
nection electric field, which was found to be proportional 
to the product of available magnetic energy density at both 
sides. The maximum of the peak reconnection electric field 
was shown to bisect the total magnetic shear angle (p. 
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The principle that determines the orientation of the x- 
line in this simple planar current sheet could potentially 
guide us to find the location of reconnection in a more real¬ 
istic magnetopause geometry. Global magnetospheric MHD 
simulations were recently performed [Komar et al.^ 2015] to 
compare these models, along with other ideas that predict 
the locations first in the global geometry with the “orienta¬ 
tion” being the resulting locus that connects these locations. 
These predictions include maximizing the total magnetic 
shear angle [Trattner et al.^ 2007], the total current density 
Alexeev et al, 1998] and the divergence of the Poynting flux 
Papadopoulos et al, 1999]. Another scenario suggested the 
X-line to be the magnetic separator simply resulting from 
the vacuum superposition of Earth’s dipolar and solar wind 
magnetic fields [Cowley, 1973; Siseoe et al, 2001; Dorelli 
et al, 2007], which was also tested [Komar et al, 2013]. 

Observationally, the location and orientation of x-lines 
have been inferred from patterns of accelerated flows [Dun¬ 
lop et al, 2011a; Phan et al, 2006; Pu et al, 2007; Seurry 
et al, 1994], and patterns of precipitating ion dispersions 
[Trattner et al, 2007] during quasi-steady reconnections. 
Statistical studies of the flux transfer events (ETEs) gener¬ 
ated by bursty reconnection [Fear et al, 2012; Dunlop et al, 
2011b; Wild et al, 2007; Kawano and Russell, 2005], and the 
global distribution of streaming energetic ion anisotropies 
[Daly et al, 1984] also provide clues. In addition, meth¬ 
ods for locally reconstructing the reconnection geometry 
[Teh and Sonnerup, 2008; Shi et al, 2005; Denton et al, 
2012] were developed, some of these methods [Shi et al, 
2005; Denton et al, 2012] could potentially take advan¬ 
tage of satellite clusters that are deployed closely, such as 
NASA’s Magnetospheric Multiscale Mission (MMS) [Bureh 
and Drake, 2009]. These in-situ observations constrain the¬ 
oretical modeling, however, it has been difficult to use them 
to distinguish between the detailed predictions of these mod¬ 
els. To accurately determine the x-line orientation in obser¬ 
vation is still a challenging and active research area in space 
study. 

In this paper, we use 3D simulations to study the ori¬ 
entation of X-lines in a given asymmetric planar geometry. 
In a similar work, Sehreier et al [2010] used 3D Hall-MHD 
simulations, where the x-line orientation is concluded to be 
consistent with that of the maximized reconnection outflow 
speed [Swisdak and Drake, 2007] or reconnection electric 
field (rate) [Cassak and Shay, 2007]. However, a study us¬ 
ing 3D fully kinetic simulations does not exist, where kinetic 
effects such as the particle streaming could be potentially 
important. 

Eurthermore, we develop a way to clearly test the orien¬ 
tation of a single x-line in 3D systems that develops in a 
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controlled fashion. Eor reconnections that develops from a 
long current sheet without a perturbation, or with a per¬ 
turbation similar to that in the GEM-challenge [Birn et al.^ 
2001], the X-line orientation may be strongly affected and 
even selected by oblique flux ropes arising from tearing in¬ 
stabilities in the linear phase (e.g., Yi-Hsin Liu et al. [2013]). 
To avoid this, we use a spatially localized perturbation to 
induce a single x-line. This localized perturbation prevents 
the linear tearing instability before the development of a 
single X-line, and does not pre-select the orientation of x- 
line. The single x-line that develops with sufficient freedom 
appears to bisect the total magnetic shear angle for larger 
mass ratios. This is consistent with that suggested in Hesse 
et al [2013]. 

The layout of this paper is the following. Section 2 de¬ 
scribes the setup of our particle-in-cell simulations. Section 
3 shows the measurements of x-line orientations in 3D simu¬ 
lations with mil Trie — 25 and mijme — 1 plasmas. Section 4 
compares the results with 2D simulations, and the mass ra¬ 
tio dependency of x-line orientation is investigated. Section 
5 contains the summary and discussions. 

2. Simulation setup 

The asymmetric configuration employed [Hesse et al, 
2013; Aunai et al.^ 2013; Pritehett, 2008] has the mag¬ 
netic profile, B = Bo(0.5 + aS')xo + Byoyo where S = 
(aitanh(z/A). This corresponds to a shear angle 0 = 180° — 
ta,n~^[(Byo/B q)/(0.5 + ai)] - ta,n~''-[(Byo/Bo)/(0.5 - ai)] 
across the sheet. The plasma has density n = no[l — a 2 (SY 
*S'^)/3] and an uniform total temperate T = 3Bq/( 87rnoa2). 
We choose ai = (a 2 = 1, then the resulting B 2 x 0 = 1.55o, 
Bixo = 0.5Bo and n 2 = no/3, ni = no. Here the subscripts 
“1” and “2” indicate the magnetosphere and the magne¬ 
tosheath sides respectively. We use a uniform guide field 
Byo = Bq then the total magnetic shear agnle 0 ~ 82.87°. 
The temperature ratio is T^/Te = 5, and the ratio of 
electron plasma to gyro-frequency is ojpelLlce = 4. Here, 
cjpe = (47rnoe^/me)^^^ and Ltce = eB^jmeC. 

In this paper, fully kinetic simulations were performed 
using the particle-in-cell code -VPIC [Bowers et al.^ 2009]. 
Densities are normalized by density no, time is normalized 
by the ion gyro-freqency Qci, velocities are normalized by 
Alfvenic speed Va = Ho/(47rnomi)^/^, and spatial scales are 
normalized by the inertia length dj = cjujpj^ where j = e,z 
for electron or ion respectively. 

Eor the rest of this paper, the x-line orientation will be 
quantified using the angle Q respect to the yo-axis. The sim¬ 
ulation box can be rotated to x = cos^boccXo YsmOboxYo and 
y = sinObox^o — cosOboxYo- In a 2D system, this machinery 
allows us study the reconnection with a pre-selected x-line 
orientation 0 = Obox- The in-plane magnetic field vanishes 
at Zn = Atanh“^{-[0.5 + (Byo/Bo)tSinObox]/oA}- 

The primary 3D run (case k in Table 1) discussed in detail 
uses mi I me — 25 and has a domain size of Lx x Ly x Lz = 
64di X OAdi X lOdi with 1024 x 512 x 256 cells. The simula¬ 
tion domain is rotated to ^ = +10°, this does not affect the 
conclusions in this paper. The boundary conditions are peri¬ 
odic both in the x- and y-directions, while in the z-direction 
are conducting for fields and reflecting for particles. We use 
150 particles per cell. The half-thickness of the initial sheet 
is A = 0.8di. In addition to 3D simulations, 2D runs with 
mijme — 1,4,25,100 and 256 are also conducted to study 
the mass ratio dependency. These runs are listed in Table 
1 . 

To study the simplest situation with a single x-line, we 
want to avoid the development of tearing instabilities be¬ 
fore a well-defined x-line forms. We use a perturbation lo¬ 
calized in the x-direction since the tearing mode is more 
stable in a short current sheet. In addition, a perturba¬ 
tion being uniform in the y-direction might pre-select the 


orientation. Therefore, we further localize the perturba¬ 
tion in the y-direction so that the single x-line can de¬ 
velop with sufficient freedom. The perturbation used in 
mijme — 25 cases is illustrated in Eig. 1. The pertur¬ 
bation has the functional form Bz oc cos[7r(z — Zp)/L^] x 
sin(27rx/Lcc)exp(-|x|/Lpi) x f(y), where f(y) = tanh[(j/+ 
Lp 2 )/Lp 3 ]-tanh[(y-Lp 2 )/Lp 3 ]. is derived using V-B = 
0 and By = 0. The peak value of the perturbation is 
5Bz — 0.05Ho. Eor runs with mass ratio mijme = 1, 4, 25 
and 256, we choose Lp\ — La,/40, LxjS, Lxj20 and Lxj20] 
Lp 2 jdi — 1, 0.5, 1 and 1; Lpzjdi = 2, 3, 2 and 2 respec¬ 
tively. To simplify the comparison, we fix Zp — — 0.5493L 
for all cases. This is the location where the in-plane mag¬ 
netic vanishes at the ^ = 0° plane. 

It may be argued that the orientation of these single 
X-lines will be constrained to meet the resonant condition 
imposed by periodic boundaries, i.e., (Ly/La,)tan(^ + Obox) 
being a rational number, which is equivalent to the safety 
factor in fusion Tokamaks (e.g., Beidler and Cassak [2011]). 
An extensive study on the effect of periodic boundary using 
mass ratio mijme = 4 is performed, but not shown here. 
We find that the x-line in a large enough simulation box 
develops into the same orientation even with different box 
aspect ratio and box orientation. We conclude that the lo¬ 
calization of X-line as described mitigates this effect from 
periodic boundaries. 

3. 3D simulation results 

With the localized perturbation, a well-defined x-line 
emerges near the center of the simulation box. Eigure 2 
shows the total current density |J| of the primary case at 
time 60/Qcz, when a single x-line with large scale outflows 
has developed. The planes are cuts at z = 0 and y = 0. In 
order to measure the orientation of the x-line, we focus on 
the X — y plane (top-view) in Eig. 3. The current density is 
shown in panel (a) and a black-dotted line of ^ = —13° is 
overlaid for comparison. To avoid a potential dependency on 
the choice of the x — y plane, the 3D iso-surface of | J| = 2 is 
plotted in Eig. 3(b), which further justifies the measurement 
of this angle. The reconnected magnetic field Bz is shown in 
Eig. 3(c). The region with Bz = 0 indicates the topological 
separator and follows the same black-dotted line. Note that 
the average magnitude of Bz is around ~ 0(0.1 Ho), as ex¬ 
pected in the nonlinear stage of reconnection. Eigure. 3(d) 
depicts the non-ideal electric field Hy, which traces the dif¬ 
fusion region of magnetic reconnection, also shows the same 
orientation. Eig. 3(b) and (d) suggest the x-line extension 

20di = lOOde, and interestingly the x-line does not ap¬ 
pear to extend much longer at later time. A similar finite 
extension was observed in symmetric reconnection simula¬ 
tions [Shay et al, 2003]. 

To evaluate the global reconnection rate, we apply the 
general magnetic reconnection theory (GMR) [Sehindler 
et al., 1988; Hesse and Sehindler, 1988; Hesse and Birn, 
1993] on this three-dimensionally localized x-line. GMR the¬ 
ory points out the importance of evaluating the integration 
of the parallel electric field Hy along magnetic field lines, 
E = f Ejjds, especially for field lines that thread the ideal 
region (Hy = 0) through the non-ideal region (Hy ^ 0) 
to the ideal region at another end. The maximum value, 
Emax = max[S(x,z)], is the global reconnection rate. This 
will be an accurate measure of reconnection rate since the 
net contribution of electrostatic component in Hy, that is 
not directly relevant to reconnection, will vanish in this in¬ 
tegration. The integration reduces the 3D system to a 2D 
map of E, as shown in Eig. 4(a). This E map in the y = 0 
plane is generated by integrating Hy along field lines for 
30di arc-length at both sides of the y = 0 plane. We can 
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then identify the location of Smacc on this 2D map and trace 
the magnetic field line from this seed point (yellow). This 
magnetic field line that carries Smacc is expected to be tan¬ 
gential to the X-line locally around the diffusion region if 
the diffusion region is quasi-2D for a reasonably long ex¬ 
tension. Eor comparison, we also trace 15 field lines seeded 
evenly along the z-direction at the same x and y coordinate 
of 'Emax- These sample field lines with positive (negative) 
Bx are colored in red (blue). 

Eigure 4(b) shows the top-view of these field lines over¬ 
laid with the iso-surface of Ey = O.OSVaBq/c (green). The 
field line with Emax (yellow) appears to pass through the 
non-ideal region and is tangential to the black-dashed line 
with orientation 0 = —13°. This orientation approximately 
bisects the total magnetic shear angle across the current 
sheet (i.e., the angle between the red and blue field lines). 
It may be argued that the field line behavior may be sensi¬ 
tive to the choice of seed points due to the chaotic nature 
of magnetic field lines [Boozer, 2012]. Hence, to get a more 
conclusive measurement, we also seed 100 points evenly dis¬ 
tributed inside a sphere of radius 0.1 de centered at the loca¬ 
tion of Emax^ These field lines traced from these seeds are 
shown in Eig. 4(c) in yellow. They align with orientation 
0 ^ —13° inside the non-ideal region (green), then separate 
quickly from each other outside the non-ideal region. The 
global reconnection rate is Emax ~ 4:.SVABodelc. Divided 
by the length of the x-line ~ lOOde, the spatially-averaged 
2D rate is roughly 0.048VaBo/c, quantitatively similar to 
the peak measured in the corresponding 2D simulation 
at this orientation (Eig. 6(c)). This further justifies that 
this 3D X-line is at its nonlinear phase. 

In Eig. 3(b), theoretical predictions of x-line orienta¬ 
tion [Sonnerup, 1974; Swisdak and Drake, 2007; Cassak and 
Shay, 2007; Sehreier et al, 2010; Birn et al, 2010; Hesse 
et al, 2013] are plotted as dashed lines using different col¬ 
ors. Note that a prediction based on the reconnection elec¬ 
tric field in Birn et al. [2010] is also presented here, where a 
more accurate energy equation is considered to improve the 
Cassak-Shay formula [Cassak and Shay, 2007]. The ratio of 
specific heats 5/3 is used. The closest prediction is the an¬ 
gle of bisection with 0 = —14.87° [Hesse et al, 2013; Moore 
et al, 2002; Borovsky, 2008; Sibeek, 2009]. Using this same 
asymmetric configuration with milme = 25, Hesse et al. 
[2013] found a relation between the peak reconnection elec¬ 
tric field and the available magnetic energy for reconnection 
Erec oc Hi^^ec^ 2 ,rec- The oricutatiou that bisects the total 
magnetic shear angle maximizes this Erec- 

While the agreement between this 3D simulation and the 
theoretical prediction is excellent, it is important to test if 
this bisection orientation is generic. We perform a similar 
3D simulation in electron-positron plasmas with mass ratio 
milnie = 1 to test the mass ratio dependency. The mea¬ 
surement of |J|, Bz and Ey displayed in Eig. 5 consistently 
suggest an angle ~ —28°, which is larger than the bisection 
angle. However, all existing analytical predictions [Son- 
nerup, 1974; Swisdak and Drake, 2007; Cassak and Shay, 
2007; Sehreier et al., 2010; Birn et al., 2010; Hesse et al., 
2013] remain the same as indicated in Eig. 5(b), since they 
do not have the mass ratio dependency. This is investigated 
further in the following section. 

4. 2D modeling and predictions 

To understand the difference with a lower mass ratio, we 
go back to 2D simulations at oblique planes (i.e., Obox ^ 0). 
Unlike 3D, the advantage of using 2D simulation is that 
we can choose the orientation of the x-line, which is out 
of the 2D plane. We study the evolution of the reconnec¬ 
tion rate Rq = {dt'ip/dt) /(VaBq) at different orientations 


in Eig. 6. Here 'ip is the difference of the flux function Ay 
between the primary x- and o-points, which are the saddle 
point and local maximum of Ay respectively. As shown in 
Eig. 6(c), with milme — 25 the orientation that maximizes 
the peak reconnection rate are consistent with the bisection 
angle with 0 = —14.87° [Hesse et al., 2013]. However, with 
a lower mass ratio the orientation that maximizes the peak 
rate shifts to a larger angle as shown in Eig. 6(a)-(b). The 
suggested angle from these 2D simulations with electron- 
positron plasmas (mijme = 1) is ^ ~ —28°, consistent with 
the orientation measured in the 3D system (Eig. 5). This 
further suggests that 2D models are sufficient to capture the 
physics that determines the x-line orientation in 3D systems, 
and this orientation maximizes the peak reconnection rate 
among these 2D oblique planes. With a larger mass ratio, 
the bisection angle persists to maximizes the peak rate in 
2D simulations, as shown in Eig. 6 (d) with mijme — 256. 
While a 3D simulation similar to that of Eig. 2 with a 
realistic mass ratio mijme — 1836 is impossible with cur¬ 
rent computational capability, the consistency between 3D 
and 2D simulations demonstrated here suggests that the x- 
line may still bisect the magnetic shear with real mass ratio 
mijme — 1836. This prediction is directly relevant to the 
reconnection events at Earth’s magnetopause. 

To explain the mass ratio dependency, we notice that sec¬ 
ondary plasmoids are generated to cause fluctuations in the 
reconnection rates shown in Eig. 6(a)-(b), even though we 
used the localized perturbation and a thicker initial sheet. 
In contrast, the induced singe x-line in plasmas of higher 
mass ratios (e.g., mijme — 25,256 in Eig. 6(c)-(d)) does 
not generate secondary plasmoids, presumably because the 
Hall effect arising from the mass ratio difference prevents the 
opened reconnection exhaust from collapsing [Shay et al., 
1999; Stanier et al., 2015], and hence makes tearing modes 
more stable This further motivates us to conjecture that 
the physics of tearing instability may play some role, that is 
more apparent with a lower mass ratio. The tearing instabil¬ 
ity is driven by the filamentation tendency of current sheet. 
In principle, the nonlinear current sheet of the single x-line 
could still be subject to the same filamentation tendency. 

To investigate this, we use a current sheet of de-scale 
thick. The half-thickness A = 1.36de is taken to be the 
mean value of the inertial lengths at both sides (i.e., de and 
VSde). This thickness mimics the current sheet scale ob¬ 
served in the nonlinear stage of reconnection. In Eig. 7 
(a), we show the tearing modes that spontaneously grow in 
this de-scale current sheets without any perturbation. At 
the time of measurement, the tearing mode amplitude is 
still small, SBzjBo ~ 0(10“^), and hence justifies its linear 
stage. The dominant tearing modes, presumably the fastest 
growing tearing mode, has a similar orientation 0 ^ —28° 
as that of the single x-line observed in Eig. 5. This sup¬ 
ports our conjecture on the role of the tearing instability. In 
addition, we do the same experiment with higher mass ra¬ 
tio mijme — 25. Interestingly, the dominant tearing modes 
manifest an angle 0 ~ —13° as shown in Eig. 7(b), that is 
also consistent with the x-line orientation in Eig 3. These 
results imply that tearing modes may have some relation 
to the peak reconnection rate measured in 2D simulations 
(Eig. 6), and hence the bisection solution in the large mass 
ratio limit. 

5. Summary and Discussion 

We demonstrate that in the large mass ratio limit the x- 
line bisects the total magnetic shear angle across the current 
sheet, at least in the 3D simulations presented here. The ori¬ 
entation can generally be predicted by scanning through a 
series of 2D simulations to find the orientation that maxi¬ 
mizes the peak reconnection electric field. This result serves 



X - 4 


YI-HSIN LIU ET AL.: X-LINE ORIENTATION 


as a practical prediction to reconnection events at Earth’s 
magnetopause. 

The fact that de-scale tearing modes share the same ori¬ 
entation as a nonlinear single x-line may have profound im¬ 
plications. The tearing instability is driven by the filamen- 
t at ion tendency of the current sheet. In principle, the non¬ 
linear current sheet of the single x-line could still be subject 
to the same tendency, and consequently develop into a state 
that is marginally stable to the tearing instability. The lin¬ 
ear tearing mode hence may provide predictions on some 
properties of the x-line, for instance, the orientation shown 
here and maybe the spatial scale of the x-line [Yi-Hsin Liu 
et a/., 2014]. To study these de-scale tearing modes in 3D 
simulations with a larger mass ratio is still computationally 
feasible, since we only need to check the linear phase and 
the tearing instability’s growth rate is large for a de-scale 
current sheet. In terms of ion gyro-frequency, the growth 
rate [Yi-Hsin Liu et al.^ 2013; Daughton et al.^ 2011] is 
^jLlci ~ (de/A)^(mi/me)(pe/de). In this Harris-type equi¬ 
librium [3 ~ 0(1), then pe/de ~ 0(1), hence the growth rate 
is proportional to the mass ratio milnie in a de-scale current 
sheet. Therefore, simulations like Eig. 7 could also serve as 
a useful indicator in predicting the x-line orientation. 

Some caveats and limitations need to be kept in mind. 
Eirst, at late times, periodic boundaries may start to affect 
and secondary flux ropes (i.e., 3D version of plasmoids) de¬ 
velop along the separatrix [Daughton et a/., 2011] of these 
single X-lines. These oblique flux ropes intertwine with each 
other and complicate the current sheet as seen in Eig. 8(d), 
where a definite measurement of the x-line orientation be¬ 
comes difficult. However, the orientation of the primary 
topological separator in Eig. 8(d) appears to remain simi¬ 
lar. Second, for multiple x-lines that develop from periodic 
tearing modes in a current sheet without a localized per¬ 
turbation, the orientation could be strongly affected by the 
nonlinear flux-ropes [Yi-Hsin Liu et a/., 2013]. Third, this 
study employs one possible asymmetric configuration where 
the stabilization by the diamagnetic drift [Swisdak et al, 
2010; Phan et al.^ 2010] is weak. The difference of plasma- 
(3 between both sides is ~ 2(^/di)tan((/)/2) ~ 2 if the 
current sheet thickness 6 ^ di is assumed. Euture work will 
explore the regime with ^ 2(^/di)tan((/)/2) to demon¬ 
strate the effects of diamagnetic drifts on the development 
of X-lines in 3D systems. 
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Notes 

1. However, the reconnection rate has the same order regardless 

the difference in mass ratio. 
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Table 1: Parameters of Runs 


Cases 

Type 

rriilme 

\/di 

Lx / di 

Ly ! di 

L z / di 

dbox 

dxline 

a 

3D 

1 

2.5 

256 

256 

32 

10° 

^ -28° 

b 

3D 

1 

1.5 

128 

128 

25.6 

0° 

-28° 

c 

2D 

1 

2.5 

256 

NA 

32 

[10°,-37°] 

dbox 

d 

tearing 

1 

1.36 

256 

256 

32 

10° 

NA 

e 

3D 

4 

1.5 

128 

128 

25.6 

0° 

-21° 

f 

2D 

4 

1.5 

128 

NA 

25.6 

[10°,-34°] 

dbox 

g 

3D 

4 

1.5 

83.2 

83.2 

25.6 

0° 

-21° 

h 

3D 

4 

1.5 

64 

64 

25.6 

0° 

-21° 

i 

3D 

4 

1.5 

64 

83.2 

25.6 

0° 

-21° 

j 

3D 

4 

1.5 

64 

83.2 

25.6 

-15° 

-21° 

k 

3D 

25 

0.8 

64 

64 

16 

10° 

-13° 

1 

3D 

25 

0.8 

64 

64 

16 

0° 

-13° 

m 

2D 

25 

0.8 

64 

NA 

16 

[10°,-33°] 

dbox 

n 

tearing 

25 

0.272 

64 

64 

16 

10° 

NA 

o 

2D 

256 

0.8 

64 

64 

16 

o 

CO 

CO 

1 

o 

O 

T-1 

dbox 
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Eigure 1: mil me — 25. The localized perturbation with amplitude 6Bz = 0.05Bo initially imposed in the simulation box 
with the y-direction aligned to ^ = 10°. This 2D plane is taken at z = Zp = —0.5493A. 



Eigure 2: Case-k with mijme = 25 at GO/r^ci- The global structure of the x-line is shown as \J\ in 2D planes at y = 0 and 
z = 0. The y-direction of the simulation box is aligned to ^ = 10°. 



Eigure 3: Case-k with mijme — 25 at 60/Dcz- In (a), the current density |J| at z = 0. In (b), the 3D iso-surface of |J| = 2 
is overlaid with theoretical predictions. In (c), the reconnected magnetic field component Bz Sil z — 0. In (d), the parallel 
(non-ideal) electric field E\\ normalized to VaBqIc is displayed at z = 0. Black dotted lines have 0 = —13°. 
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Eigure 4: Case-k with rmlme — 25 at 60/Qci- In (a), the 2D map of quasi-potential S = /Eyds at y = 0. The field line 
with Smacc is coloied in yellow and sample field lines distributed vertically are colored in red {Bx > 0) and blue (Bx < 0). 
Panel (b) is the top view of (a) with the iso-surface of E\\ = O.OSVaBq/c overlaid in green. Panel (c) has a similar format 
of (b), but depicts 100 magnetic field lines (yellow) traced from seeds that are evenly distributed inside a sphere centered 
at the location of Smax with radius O.lde. 



Eigure 5: Case-a with milrrie = 1 at 300/Qcz- In (a), the current density \J\ at z = —0.4de. In (b), the 3D iso-surface of 
I J| = 0.62 is overlaid with theoretical predictions. In (c), the reconnected magnetic field component at z = —0.4de. In 
(d), the iso-surface of parallel (non-ideal) electric field = O.OIGVaBq/c. Black dotted lines have 0 = —28°. 
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Eigure 6: Reconnection rate Ro at sample orientations with mass ratios (a) mil me — 1, (b) mijme — 4, (c) mijme — 25 
and (d) mijme — 256. 


(a) mil me = 1 (b) mijme = 25 



Eigure 7: Tearing modes in a electron-scale current sheet of half-initial thickness A = 1.36de. In (a), case-d with mijme — 1 
at time 51.25/Qci- Rz is time-averaged for 12.5/Qci- The black dotted line has Q — —28°. In (b), case-n with mijme — 25 
at time 3.5/Qcz- Rz is time-averaged for 2lQ.ci- The black dotted line has 0 — —13°. 
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Eigure 8: Case-k with rriilme = 25. The iso-surface of the reconnected magnetic field component BzjBo = 0.08 is shown 
for time (a) t = 40/Qcz, (b) t = 60/Qcz, (c) t = 80/Qcz and (d) t = 100/Qcz- The black dotted lines have 0 = —13°. 































